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The Fock space of a system of indistinguishable particles is isomorphic (in a non-unique way) to 
the state-space of a composite i.e., many-modes, quantum system. One can then discuss quantum 
entanglement for fermionic as well as bosonic systems. We exemplify the use of this notion -central 
in quantum information - by studying some e.g., Hubbard, lattice fermionic models relevant to 
condensed matter physics. 
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A big deal of attention have been recently devoted to 
the notion of quantum entanglement [jjj. As a matter 
of fact this fashion is mostly due to the vital role that 
such notion is generally believed to play in Quantum In- 
formation Processing (QIP) protocols S. The definition 
of entanglement relies on the tensor-product structure of 
the state-space of a composite quantum system ||. 

On the other hand such a tensor product structure is 
not present in a large class of systems of major physical 
interest: ensembles of indistinguishable particles. Indeed 
in this case it is known - basically since the birth of quan- 
tum theory- that the state-space associated with N sub- 
systems is constrained to be subspace of the A-fold tensor 
product. Depending on the bosonic or fermionic nature 
of the subsystems one has to select the totally symmetric 
or antisymmetric subspace. This is an instance of a su- 
perselection rule || i.e., a fundamental limitation on the 
possibility of preparing a given state. 

We see therefore that the existence of quantum statis- 
tics H makes the notion of entanglement problematic 
for systems of made of indistinguishable subsystems e.g., 
particles. Notice that, even if one allows for general 
parastatistics ||, most of the product states ofljji) do 
not belong to the physical state-space at all in that they 
do not have a proper transformation under permutations 
of Sn i.e., they do not belong to an ^Ar-irrep 0]. 

Very recently some authors addressed the issue of en- 
tanglement (or more generally quantum correlations) in 
system of two fermions Q and bosons . Their approach 
appears to be the natural generalization of the one used 
for the distinguishable particles. 

In this paper we shall tackle the problem of the rela- 
tion between entanglement and quantum statistics from 
a rather different perspective based on entanglement rel- 
ativity as discussed in Ref. Q We shall mostly focus on 
the fermionic case j^]. In particular we shall analyze the 
local i.e., on-site, entanglement associated with simple 
fermionic models on a lattice. 

Fermions and Qubits. Let us start be recalling ba- 
sic kinematical facts about many fermion systems. Let 
h L := span{|V>z)}ieN L (N L := {1, ...,£}) be an L- 
dimensional single particle state space. The labels in 
Ni will be referred to as sites and the associated single- 
particle wavefunctions will be thought of as describing a 
(spatially) localized state. Accordingly the set of the i's 
will be referred to as the lattice. 



Let Vl (T-'l) denote the whole (with N elements) fam- 
ily of subsets of Nl. For any A := {ji, . . . , j^} G V^ we 
define the anti-symmetrized state-vector 
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The |A)'s are an orthonormal set. The state-space 
T-Ll{N) associated with TV (spinless) fermions with single- 
particle wavefunctions belonging to is given by the 
totally anti-symmetric subspace of hf N i.e., Hl(N) := 
span{|A) / A S V£}. The fermion number ranges from 
to L, the total Fock space is obtained as a direct sum of 
the fixed number subspaces i.e., Hl = ®n = qHl{N) — 
span{|A) I A € Vl}- From the well-known relation 

dimH L = J2N=o dimH Li N ) = ELo ( n 

lows that the fermionic Fock space is isomorphic to a L- 
qubits space, each qubit being associated with a site ]To| |. 
The latter isomorphism is realized by the the mapping 



= 2 L ,it fol- 
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where xa-^l — * {0, 1} is the characteristic function of 
A. Clearly A(|A)) is nothing but a A-qubit basis state 
having in the j-th site a one (zero) if j € A (j ^ A) . The 
0-particle state |0) is mapped by A onto |0) := \0)® L ; 
thus latter vector is referred to as the vacuum. 

In our considerations, once TIl is endowed by A with 
a multi-partite structure, tensor products of individual 
single-particle spaces are not relevant anymore. To exem- 
plify this point let us consider the case L = 3. It is not dif- 
ficult to see that all the states in 7i 3 (2), seen as elements 
of hf 2 , have the same entanglement. Indeed all of them 
can be written as |o) ® |6) — |6) ® |o), for suitable |o) and 
| b) |pH . On the other hand both the "separable" state 
|1)®|1)®|0) and the "entangled" (|0)®|1)-|1)®|0))®|1) 
belongs to A (7^3(2)). This kind of puzzle is solved by ob- 
serving that the entanglement of, say \a) <g> \ b) — \b) ® \a), 
is not physical. Indeed in the involved subsystems i.e., 
individual "labelled" particles, due to the very notion of 
indistinguishability, are physically not accessible. 

This situation is just an illustration of the relativ- 
ity of the notion of entanglement ||. The latter cru- 
cially depends on the choice of a particular partition 
into physical subsystems. In this case "good" subsys- 
tems are associated with the set of single particle modes 
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(labelled by I £ whose occupation numbers are phys- 
ical observables and not with the particles themselves. 
From this perspective one can have entanglement with- 
out entanglement. For instance a one-particle state e.g., 
|0) <g> |1) + |1) <£> |0), can be - with respect to the parti- 
tion into mode subsystems- entangled. It is important 
to stress that such kind of one-particle entanglement (or 
correlation), despite its paradoxical nature, has been re- 
cently proven to allow for quantum teleportation [ ^2[ ; 
therefore it has to be regarded as a genuine resource for 
QIP. 

The Fock space TL l , since it allows for a varying parti- 
cle occupation, does not correspond generally to the state 
space of a physical system. For charged fcrmions coherent 
superpositions of vectors belonging to different particle 
number sectors are forbidden due to charge superselec- 
tion rule || . In this sense our qubits are unphysical. Only 
qubit states in the A(Hl(N)) are associated with (N- 
particle) physical states. Accordingly not all the elements 
of End(H.L) correspond to physical observables: the lat- 
ter span the subalgebra T of number conserving opera- 
tors i.e., T := {X I [X, N] = 0} = ® N End(H L {N)) 

Local Entanglement .Let |\&) S Hl(N) be the associ- 
ated j-th local density matrix is given by pj := Tvj |\&) (W| 
where Tr, denotes the trace over all but the j-th sites. 
For any j £ N^, one obtains a bipartition of Hl i.e., 
C 2 (g> (<D 2 )® therefore the entropy S (von Neumann 
as well as linear) of pj is a measure of the entanglement 
of the j-th site with the remaining N — 1 ones. 

Local entanglement is relative to the decomposition 
into subsystems i.e., sites, defined by mapping (||) ||. 
One could consider different isomorphisms giving rise to 
inequivalent partitions in to "local" subsystems. This 
fact can be clearly seen by introducing creation and an- 
nihilation operators {cj}f =1 C End(Wt), ( [a, b]± = 
ab ± ba) which satisfy canonical (anti) commutation re- 
lations for (fermions) bosons 

h, Cj}± = 0, [a, c]]± = Cj |0) = (j e N L ). (3) 

Of course Hl = span{n^ = i(cj)" J |0) / n x , . . . , n L £ 
Nqo}. If U is a L x L unitary matrix then it is well-known 
that the following (Bogoliubov) transformation 

L 

Ci -> h := Uij c j: (i £ N L ) (4) 

2=1 

maps fermions ( bosons) onto fermions (boson) giv- 
ing rise to an automorphism of the observable alge- 
bra. Accordingly new occupation-number representa- 
tions Aj/: nf=i (cl) ni |0) ^ ®f =1 \n t ) (n z = 0,1) are de- 
fined. Clearly entanglement is strongly relative to the de- 
compositions associated with different Ay's. Notice that 
even though automorphsims (^) have a single-particle ori- 
gin, they define non-local transformations of the Fock 
space onto itself. Indeed a mapping W £ U((Hl) is lo- 
cal with respect the subsystem decomposition associated 
with A v iff A v o W o Ajj 1 £ nf =1 U{€ 2 )i 



As a particular, though quite relevant, case one 
can consider the Fourier transformation i.e., U^j ■= 
e ikj , k := 2n(l - 1)/L,(l £ N L ). The wave-vectors k 
label the so-called reciprocal lattice (Ajj is denoted by 
A*) and represent physically modes delocalized over the 
spatial lattice. It is obvious that states that are entangled 
(non-entangled) with respect A can be non non-entangled 
(entangled) with respect A*. 

The situation we shall investigate in this paper is the 
following. Suppose H £ End(7ii) is a non degenerate 
(gran-canonical) hamiltonian and H — ^2 m e m \^m)(^m\ 
its spectral decomposition. If pj 1 denotes the j-th local 
density matrix associated with the energy eigenstate s m , 
one can compute the quantity 

m— 1 j— 1 

where Z is the (gran-canonical) partition function i.e., 

Z := Y^m=x e ^ em '• © ^ s ^ ne thermal expectation 
value of the local entanglement averaged over the whole 
lattice . In particular we will be interested in the limit 
f3 i * oo i.e., local entanglement So in the ground state. 

When energy spectrum shows degeneracies Eq (||) is 
no longer well-defined. We assume that there is a "nat- 
ural" (see examples below) way to select a a complete 
set of commuting observables containing H, whose joint 
eigenvectors provide the e m 's to be used in (||). 

To begin with we observe that 

Pi = |1)(1| (¥|ni|¥) + |0> <0| - n,-|*) (6) 

where ni :— cj Cj = ® lj is the local occupation 

number projector. Indeed: (l]pj|l) = Tr(|0)(0| pi) — 
Tr(rii pi) = (*|nj|*) in the same way one obtains 
the other diagonal element of pi. Moreover (0|pi|l) = 
Tr(|l)(0|pi) = Tv{c\ Pi ) = (*|4|*) = 0, last equality is 
due to the fact that |*) is a particle number eigenstate 
i.e., an eigenstate of the operator N :— nj. 

Itinerant Fermions. We now consider free (spinless) 
fermions hopping in the lattice. The Hamiltonian is given 
by 

L-l 

Hpree = ~t ^ (c} +1 Cj + h.C.) - pN (7) 

2=1 

Introducing the Fourier fermionic operators Ck ■= 
l/V^S_y=i e%k ^ c ji if is a text-book exercise to prove 
that (0) has eigenstates given by the Af-particle vec- 
tors |k) := nm=i4 ro l°)' ( k : = (h,-..,k N ) £ R N with 
eigenvalues := — 2i^m=i cos(k m ) — pN. 

The local density matrix is easily obtained by us- 
ing Eq. (|(|) and the translational properties of the 
|k)'s. If T denotes the natural representation in Hl 
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of the cyclic permutations i > i -+ 
translation operator, one has T|k) 



1 mod L i.e., the 



«=1' 



Therefore (k|n 3 |k) = 1/L Ei=i( k Kl k ) = = : 
n. Whereby E = 1/Z£^ =0 S(N/L)e^ N Z N (f3) = 
l/ZTv(S(N /L)e- f3Hp ^), in which S(n) = -n Inn - 



(1 - n) ln(l - n) and Zjv(/3) 



Tr 



H L (N)Z 



the (iV-particle) canonical partition function. 

The fraction p(N) := e®^ Zn jZ gives of course the 
probability of having any A^-particle configuration. In 
the thermodynamical limit (iV, L > oo, N/L =const) 
p(N) becomes strongly peaked around the expectation 
value Nq of N. In this case local entanglement is simply 
given by the Shannon function E ~ S(no), it readily dis- 
plays an intuitive feature: local entanglement vanishes 
for the empty (fully filled) lattice being the unique asso- 
ciated state given the product |0) moreover E is 
maximal at half-filling i.e., no = 1/2. Notice that for the 
states |k) entanglement associated with the A* partition 
is obviously zero. 

Spin 1/2 Fermions. Here we consider the lattice mod- 
els of spin 1/2 fermion model. We have then to introduce 
an extra dicothomic variable a —\, | to label the single- 
particle state- vectors. As usual fermionic operators cor- 
responding to different cr's always anti-commute. In this 
case it is convenient to consider the 2 2 L -dimensional Fock 
space as isomorphic tch L-fold tensor power of four- 
dimensional space i.e, TLf — (<D 4 )® L . The local state 
space is spanned by the vacuum |0) and the vectors 

I T)j := c] T |0), | l)j := cJjO), | Uh ■= 4^0). (8) 

The pj = Tr 3 -|*)(*| is now a 4 x 4 matrix. If the N- 
particle state |^) is a) translational invariant, b) eigen- 



state of S z 



Uji), is easy to see that 



pj = l/idiag(l — N-[ — Ni — N t , N h N h N t ) where 



N * : = E^i^K^l - »j-«r)l*>. (°" =T4) is the num- 
ber of lattice sites singly occupied by a tr fermion and 

Ni := Ylj=ii i &\ n 3f n 'jl\'®) i s the number of doubly occu- 
pied sites. We see that local entanglement, in state \$f) 
is a function just of the occupation numbers N a , (a =f 
, 1,1), in particular it follows that Eq. (Q) can be effec- 
tively computed for Hamiltonians commuting with the 
iVs i.e., E = J2 {Na } e^ N S({N Q }) Z({N a })/Z. An in- 
stance of this case is illustrated in the following. 

Super symmetric Dimer. We consider here a two-site 
i.e, a dimer, version of the so called supersymmetric EKS 
model fT4[ ]. For zero chemical potential i.e., half- filling 
the EKS Hamiltonian acts on the basis states as follows 



H\a) 



(-1) 



M t/3| 



(9) 



where \a\ is the parity of the single-site state \a) i.e., 
| t | = | i | = 1, 10] = | U I = 0. Since H is just a 
graded permutator the relations [H, N a ] = hold true. 
The state-space split according the N a 's configurations 
= ®{N a }T~t({N a }) and the Hamiltonian can be diag- 
onalizcd within each sector. Notice that Eq. M) is also 



invariant under a global particle-hole transformation i.e., 

\v)++\-<t),\0)»\U)- 

It is straightforward to check that H admits four sin- 
glets non-entangled (the configurations (0,0,0), (2,0,0) 
along with their particle-hole conjugates) and six dou- 
blets ( (1,0,0), (0,1,0) and conjugates and the self- 
conjugated (1, 1, 0), (0, 0, 1)) with entanglement ln2. 
Moreover since H 2 = 1 one gets an energy spectrum 
given by {—1, 1}, being both the eigenvalues 8-fold de- 
generated. Therefore E — 12 In 2 cosh (3/16 cosh/? = 
3/4 In 2 : the local-entanglement (at half filling) is tem- 
perature independent. 

This very simple result is due to the large symmetry 
group of the Hamiltonian (||). A more interesting case 
is obtained introducing a model in which a free param- 
eter controls the competition between the localized and 
itinerant nature of the lattice fermions. 

Hubbard Dimer. If Hp ree simply denotes Eq. (^) with 
extra spin index then Hubbard Hamiltonian reads (llfj) 



H 



Hubb — 



m 



.7=1 



(10) 



The new local terms added account for the on-site in- 
teraction e.g., Coulomb repulsion, experienced by pairs 
of (opposite) spin fermions sitting on the same lattice 
site. By introducing the total fixed spin number opera- 
tors N a := Ylj=i n j<? ( a = T> I) i s eas Y to check that both 
of them commute with the Hubbard Hamiltonian (|ic|). 
This implies that H^ubb can be separately diagonalized 
in each joint eigenspace 7i(N^,Ni) of the A^'s. In the 
the dimer case i.e, the dimer, one finds dim Ti(N^, N\) — 

Y[ a J ; then at most (for Nf = N± = 1) one has to 

solve a four-dimensional diagonalization problem. 



So 




FIG. 1. Entanglement of the Hubbard dimer ground state 
as a function of {7/4 1 for decomposition associated with real 
and reciprocal lattice. 

The (unnormalized) ground state, for the repulsive 
case U > 0, is given by —Go |0) where 



Go 



C 1T C U C 2T c 2i a +(^/^) ( C lt C 2 



-2| L 21 
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C U C 2f) 



(11) 
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where a±(x) := x ± \/T+?, and the associated eigen- 
value is given by Eq = — 2icv_. The entanglement of 
the state ( pi] ) is easily studied as a function U/At. Us- 
ing liner entropy as an entanglement measure one finds 
S (U/t) = 1-Trp 2 , = l-l/2(a| + l) (c^ + l)" 2 . Local 
entanglement is monotonically decreasing as a function 
of U /At. (Fig. 1). In particular one obtains the free limit 
So(0) = 3/4 and the strong coupling limit S'o(oo) = 1/2 
that correspond to ground states given by uniform super- 
positions of respectively four and two states [see Eq. (|ll[ . 
Of course the physical interpretation is quite simple: the 
higher the on-site repulsion U the more local charge fluc- 
tuations are suppressed and the smaller the number of 
available states. Eventually for infinite repulsion doubly 
occupied sites get decoupled and only spin fluctuations 
survive. In this regime Hubbard model is known to be 
equivalent to an anti-ferromagnetic Heisenberg model for 
spin 1/2 Ground state as well as thermal entan- 

glement for this (and related) models have been quite 
recently studied Jl7| ] . 

It is instructive to write the dimer ground state cre- 
ator ([ll]) in terms of the Fourier operators c„ := 
l/y/2 {da ± c 2(T ) (o- =t, I), from Eq @ one finds 



Go (a) 



fe=± 



(1 + ka + ) 



J*0t J*0t 
T 1 ' 



(12) 



With respect this reciprocal decomposition local entan- 
glement is an increasing function of U /At. From the free 
case (a+ = 1), that is unentangled up to strong coupling 
(a + = oo) which gives 5q(oo) = 1/2 (see Fig. 1). 

The example of the Hubbard dimer shows that -not 
surprisingly - entanglement is well-suited to analyze the 
interplay between itinerant and localized features of Hub- 
bard Model (|l0|): hopping term t (repulsion U ) term is 
responsible for entanglement in the real (reciprocal) lat- 
tice decomposition. 

Conclusions. In this paper we discussed some issues 
related to entanglement in system of indistinguishable 
particle. For these systems quantum statistics applies 
and therefore their state-space is not naturally endowed 
with a tensor product structure. 

Nevertheless mappings between their Fock spaces and 
multi-partite state-spaces can be established (the well- 
known occupation number representation) and then the 
usual definition of entanglement can be applied. For sys- 
tems with L single particle states available the set of 
possible inequivalent decompositions into L subsystems 
(modes) is parametrized by the group U(L) of Bogoli- 
ubov transformations. 

We focused on simple e.g., Hubbard, models of 
fcrmions on lattice studying how, as function of the model 
parameters, local entanglement varies both with respect 
to real and reciprocal lattice decomposition. Results sug- 
gest that this notion of entanglement is well suited to 
describe interplay between localization and itinerancy in 
these systems. 



We believe that the approach pursued in this paper 
- besides to establish a connection between the field of 
Quantum Information Processing and condensed matter 
physics - can provide novel physical insights in the study 
of interacting ensembles of indistinguishable particles. 
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Ref. || for drawing my attention on their work and for 
stimulating correspondence 
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